Brane-Induced Gravity's Shocks 
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We construct exact gravitational field solutions for a relativistic particle localized on a tensional 
brane in brane-induced gravity. They are a generalization of gravitational shock waves in 4D de 
Sitter space. We provide the metrics for both the normal branch and the self-inflating branch DGP 
braneworlds, and compare them to the 4D Einstein gravity solution and to the case when gravity 
resides only in the 5D bulk, without any brane-localized curvature terms. At short distances the 
wave profile looks the same as in four dimensions. The corrections appear only far from the source, 
where they differ from the long distance corrections in 4D de Sitter space. We also discover a new 
non-perturbative channel for energy emission into the bulk from the self-inflating branch, when 
gravity is modified at the de Sitter radius. 



PACS numbers: ll.25.-w, 11.25.Mj, 98.80.Cq, 98.80.Qc 

In brane-induced gravity gravitational force be- 
tween masses on the brane should arise predominantly 
from the exchange of an unstable bulk graviton reso- 
nance, which falls apart at very large scales, instead 
of a stable zero-mode graviton governed by the usual 
four-dimensional (4D) General Relativity (GR). Far from 
sources the law of gravity is modified, and changes from 
4D to a fully higher-dimensional one. Subtleties of gravi- 
tational perturbation theory in DGP have been explored 
in Qi Qi IS with the particular attention on the cou- 
plings of the longitudinal gravitons. It has been sug- 
gested that the brane extrinsic curvature may play the 
role of a coupling controller: the extrinsic curvature, re- 
sponding to a source mass on the brane, may help to tame 
perturbation theory 0, . Exploring this interplay be- 
tween the curvature and the graviton couplings on exact 
solutions of DGP field equations could clarify the status 
of effective 4D theory. The DGP equations are, however, 
notoriously difficult to solve for localized sources . 

In this Letter we present the first example of an ex- 
act solution for a localized particle in DGP models. We 
derive the gravitational field for a relativistic particle 
on a tensional brane in 5D. It generalizes 4D gravita- 
tional shock waves 0,11,01 in de Sitter space [lOllllLll^ . 
We present the closed form metrics for both the normal 
branch and the self-inflating branch DGP braneworlds 
[l3|. The scalar graviton is absent in these metrics. 
One expects this in weakly-coupled perturbation the- 
ory, where the scalar graviton decouples from ultra- 
relativistic particles since to leading order its source van- 
ishes, ^ = 0, and the corrections are suppressed by the 
source mass to momentum ratio m/p <^ 1 (e.g., hence the 
gravitational shock wave of 4D GR is also a solution in 
Brans-Dicke theory). Our results show that this persists 
in DGP, and that the dangerous strongly coupled mode 
found in perturbation theory does not destroy the solu- 
tions at large momenta after all. Further, we find that 
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as long as the brane-localized terms are present, at short 
distances the wave profile behaves exactly like in 4D GR 
0, IB| ■ The deviations appear only far from the source, 
where they differ from the 4D GR corrections in de Sitter 
space. We also suggest a new way for checking graviton 
effective field theory with our shock waves. 

DGP braneworlds fUl are given by a bulk action with 
metric kinetic terms in both the bulk and on the brane, 
with couplings set by bulk and brane Planck scales, K4 = 
1/A/| and — 1/M|. In 5D, the field equations are 



where 



5 B 



are bulk and 3-brane-localizcd Einstein 



tensors; T^^ is the brane stress-energy; {A, B} and {/i, v} 
are bulk and 3-brane world- volume indices, respectively. 
We use Gaussian normal coordinates for the 5D metric, 
ds| — dw^ + g ^j_v{'w , x)dx^ dx^ , orbifolding by a Z2 sym- 
metry w —w around the brane at w — 0. For a brane 
with tension A ^ and no matter, T^^^, = —\5^,j. The 
solution is a de Sitter 3-brane in a flat bulk , 



dsl = (1 - eH\w\fdsl^s + dw"^ 



(2) 



Apart from e = ±1 in l^), this is identical to the 5D ver- 
sion [31 of the inflating VIS domain wall 15]. We shall 
use ds|^5 = —{l — H'^r'^)dt^ + jY^jj2prj+T'^dVL2, the static 
de Sitter metric instead of the spatially flat one employed 
in p^. The space-time of (0) is a 4D de Sitter hyper- 
boloid in a 5D Minkowski bulk. Here e arises because 
we can retain either the interior e = +1 or the exterior 
e = — 1 of the hyperboloid after orbifolding, thanks to 
the brane curvature. The junction conditions relate the 
4D curvature and the brane tension as 1131 
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(3) 
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The solutions with H > differ for each of e = ±1 {H < 
cases are their PT transforms). On the normal branch, 
e = +1, the solution reduces to the 4D Friedman equation 
in the hmit M5 < M4, ~ A/Af|. In the bulk, we 
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keep the interior of the hyperboloid. It has finite volume 
and so the perturbative 4D graviton exists. The brane 
curvature terms ~ / d'^Xyi^R4/2 suppress the cou- 
pUngs of the > graviton KK modes, producing the 
4D effective theory when M5 <C M4. On the self- inflating 
branch with e = -1, at low tensions A < 12M|/M|, the 
eq. ® yields H ^ 2M|/M|. In the bulk we keep the 
exterior of the hyperboloid which has infinite bulk vol- 
ume and so there is no perturbative 4D graviton. The 
effective theory arises only from the exchange of the bulk 
resonance. 

Suppose now there is also a single relativistic parti- 
cle, say a photon, with a momentum p on the brane, 
moving along some null geodesic of the brane-induced 
metric. Its stress-energy tensor sources an additional 
gravitational field. To find it, one could boost the lin- 
earized gravitational field for a massive particle by an 
infinite amount, simultaneously taking the limit m —^ 
such that m cosh 7 = p stays finite ,7, 10, 11, 16]. View- 
ing the metric as an expansion in the powers m, in the 
relativistic limit the terms of order higher than m van- 
ish because there is only one factor of the boost param- 
eter cosh 7, and the linearized solution becomes exact! 
Alternatively, one can use a variant of a very elegant 
cut-and-paste technique developed by Dray and 't Hooft 
js^ for flat 4D backgrounds, and later applied to gen- 
eral 4D GR backgrounds by Sfetsos JJJ. We will fol- 
low this route here, by extending the technique to DGP. 
First, we cast Q in suitable null coordinates, defining 



" = Ti^Jl^^Mm and V = iy'i^exp(-i7t). 

We also change w to \z\ — —-^ ln(l — €H\w\) (note that 
6(w) = S{z)), to simplify evaluating the curvature ten- 
sors with a conformal map as in |l7| . The metric P|) 
becomes 



dsg = e 



eH\z\(_ 

1(1 



Adudv 



■ + ( 



,1 - H'^uv^^dn 



(4) 

Let the photon move along the w-axis on the brane, i.e. 
along the null geodesic u — 0. The observer on the North 
Pole of the brane de Sitter space sees a photon stream- 
ing away along the past horizon, starting from r = at 
infinite past. Using the trick of Dray and 't Hooft to 
construct the photon's gravitational field, we introduce a 
jump in the v coordinate at u = |^: let w — > w -I- Q{u)f 
and dv ^ dv + Q{u)df where / is the wave profile which 
depends only on the spatial transverse coordinates (in 
our case, the angles on the 2-sphere and z), and Q{u) is 
just the usual step function. Changing the coordinates 
to V = v + Q(u)f yields dv — > dv — 5{u)fdu. Substituting 
v^dv — > v,dv — 5{u)fdu in (4) and dropping the carets 
gives 
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Adudv 'iS{u)fdu^ 
[l + W^uvf ^ (1 + H'^uvf 



+ 



.f_l+dz'j, (5) 



1 + H^uv' 

Substituting (|SJ| into including the photon source in 



T'^, = ~X5^', + 2{p/ ^)guJ{e)5{4>)5{u)5'^5l, and eval- 
uating the curvature, bearing in mind that 5{u) and its 
derivatives are distributions (so u5{u] = 0, u^5'^{u) — 
and f{u)5'{u) = -f'{u)5{u) H [H), we find the inde- 
pendent field equations 18]. One is just eq. The 
other comes from the G5 component of |^ and yields 
a linear field equation for the wave profile: 



Ml 



2p 



{A2f + 2f)S{z) 



^^,S{n)Siz) . (6) 



Here A2 and S{Q) are the Laplacian and the (5-function 
on a 2-sphere, peaked at 6* = 0. In the limit M5 — > the 
bulk derivatives disappear, one can factor out 6(z) and 
recover the 4D de Sitter equation of 10, 12]. 

It is simpler to solve eq. for two sources with the 
same momentum p, runnin g in the opposite directions 
in the static patch [Tol lllj. In this way we avoid 
some unphysical divergences when the last term in ((7|) 
vanishes for / = 1, which cancel out anyway by sym- 
metry. The two-source solution correctly represents the 
limit of infinite boost of the Schwarzschild-de Sitter ge- 
ometry JT]. We add an extra term, j^S{il,')S{z), 

on the RHS of ® where 5{n') is peaked at 6* = tt. To 
return to a single source, we can take the solution and 
multiply it by 8(7r/2 — d) as in So one particle 

moves along 6 — and the other along 9 = tt. With this 
choice of trajectories we use the addition theorem for 
spherical harmonics to replace them with Legendre poly- 
nomials Pi{co8 9) in the expansion for the wave profile: 



f ^ j:Zoifi^\^)Piicos9) + >{z)Pi{~ cose)). Here 
fi^\z) are the bulk wave functions; //^'' is sourced by 
the photon at 9 = and by the photon at = tt. By 
orthogonality and completeness of Legendre polynomials, 
the modes f^^\z) obey the same differential equation, 

dlfi - 3eHd\,\fi + H\2 - l{l + l))fi - 



MiH^ f{2l + l)p 
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{2-lil + l))fi)6{z).{7) 



We interpret i5-function on the RHS by pillbox integra- 
tion as a matching condition for the first derivatives of 
fi^^ on the brane. The remaining boundary conditions 
come from requiring orbifold symmetry fi{—z) = fi{z) 
and square integrability in the bulk, so that the solu- 
tions are localized on the brane. Since both //^'' solve 

the same boundary value problem, fj^^"^ = fi — fi- 
By O; /('s are simple exponentials, and to be local- 
ized on the brane (i.e. square-integrable) they must be 
- e-[2'+(i-3e)/2]H|.|^ Because Pi{-x) = (-l)'P/(a;), the 
solution will be an expansion only in even-indexed poly- 
nomials P2i{cos9), like in 4D E3, E3- Solving for 
the bulk wave functions /;'s [13 yields the solution for 
two relativistic photons on the brane moving in opposite 
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directions: 



Al + l 
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1 



(i-^)e 



(i+e)g 



xe-[2'+(i-3^)/2l^l-lp2/(cos0), (8) 

where g = 2M|/(M|7?) = l/(-ffrc) (see pj). The solu- 
tion can be checked by direct substitution into 0- 
When g = 0, at z = this reproduces the 4D series solu- 
tion of |lO| for both e = ±1. We can bring the series © 
to a more compact form. Factorizing the coefficients of 
the expansion, and defining r = e^^'^' and x = cosO, 



(3 + eg)7rM2 ^ 2? - f 



[3+(l 



2i+(l-3e)/2 



P2/(x) . (9) 



Note that while for e = 1 the series Q is finite, for 
e = — f, the self-inflating branch solutions, it displays a 
spectacular behavior as g — > 1, when the brane tension 
vanishes, where the I = term of the first sum has a 
pole (g = 3 is regular, as seen from ©). A closer look 
jl8l | shows that in this limit the only finite I = solu- 
tion of {Tj) is the delocalized bulk mode, not present in 
(O. For g = 1 gravity is modified at the scale equal to 
the cosmological horizon on the brane. Brane and bulk 
begin to "resonate" : in a time-dependent problem, as 
g — > f slowly, the I — mode would grow larger, and 
at some point it may begin to produce delocalized gravi- 
tons, taking energy from the self-inflating brane into the 
bulk. This indicates an onset of a dramatic new instabil- 
ity, which warrants closer investigation. 

Now, recall the generating function for Legendre poly- 
nomials, (f — 2a;T -I- r^)~-^/^ = J2ilo ^^i^)''''' ■ Using 



lim,^o /; dM'+^ = - hm.^o /' ddd'+^ to 



1+13+1 



regu- 



late divergences when I + f3 < —1, integrate over t; after 
straightforward manipulations eq. lO becomes 



l-3e 
T 2 



2nMi 1 - l^g 2(3 + eg)7rM| 



X / dd{ , + ^ 

/o V Vl - 2x1? + Vl + 2x1? + 



(10) 



■[3-(l-6)g]^^^^ 



[3+(l + e)g]i9' 



Since H1Q(I is a Green's function in the transverse direc- 
tions, it has short-distance singularities at x = ±1, but is 
well-behaved elsewhere. We will not attempt to evaluate 
it for the general case. Here we only look at the formula 
for / on the brane z — (i.e. r = 1) at transverse dis- 
tances TZ well inside the cosmological horizon, TZ ^ H^^. 
Note that if we rewrite eq. © as (l-f eg)i/^ = A/(3M|) 
we find that on the normal branch g can vary between 
zero (4D) and infinity (5D), with the "resonance" men- 
tioned above at g = 1. On the self- inflating branch, g < 1 



or Tc > i/H, as long as A > 0. In the limit g ^ 1 on the 
self- inflating brane the tension vanishes. 

When g = 0, (|10|) is identical to the 4D case on both 
branches, as remarked above. Indeed, the integrals give 

The metric transverse to the null particle on the brane is 
ds2\z=u=o = dQ.2/H^ and so the proper "radial" distance 
is measured by the polar angle 6. For small angles, TZ ~ 
B/H, X = 1 — H'^TZ^ /2, and l(TT)) reduces precisely to the 
fiat 4D solution 0, Q: up to OiTZ'^/H-^) corrections, 
and with the sign conventions of fldi , 



TZ 



(12) 



The constant term in this equation can always be recov- 
ered in the 4D flat case by a diffeomorphism. The inte- 
grals for 0{g) corrections contain terms like J dC,hi{l + 
hC,) /C, and cannot be written in closed form ■ However, 
at short distances the singular ©(g) powers in the inte- 
grand precisely cancel out. In the remainder, since the 
leading order singularity in 1)10(1 is logarithmic around 
g = 0, the additional logarithms soften the corrections 
further, and render them finite as 7?. — > 0. Hence at 
short distances the solution looks exactly the same as in 
4D. 

In fact this persists for any finite value of g. If we 
rewrite the integrals (|10|l using the Euler substitutions 
C = Vl T 2x1? -I- -(- '§ and consider the limit 1 — x ^ 1 , 
we extract the leading singularity. It comes from the term 
P /i+V2(i-) ^ P In ^ + finite terms 

for all values of g and e. Substituting x = 1 — H^TZ^ /2, 
we recover the leading logarithm in ((12|l . The subleading 
corrections in general differ from the terms in the expan- 
sion of the 4D de Sitter solution 111() around the short 
distance limit in flat space (|12|) . In that case one finds 
corrections to come as even powers of HTZ. In DGP, how- 
ever, the corrections in (|13|l will come as odd powers of 
HTZ as well, starting with the linear term, signalling the 
hidden fifth dimension. Thus in general we will have 

m) - ^((l + a.i/^7^^ + ...)ln(^) 



+ const + biHTZ + h2H^TZ^ + ...j ,(13) 

where the coefficients , bk are numbers that can be com- 
puted explicitly for given values of g and e. 

To see how 5D gravity reemerges, on the normal branch 
we can take the limit g oo. Then the background re- 
duces to the inflating brane in 5D Minkowski bulk [T^ . 
Using to find the contribution of the second sum, 

and de- 



which is oc Vl, nP2i(x) = , H , 

ducing the contribution of the first sum from the integral 
(|10|l . we can write down the solution f^o in closed form 
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[l8l |. At short transverse distances x = 1 — H'^'R? /2 the 
leading order behavior of f^o^ using M|g — 2M^/H and 
ignoring a constant and higher powers of TZ, is 

Notice the relative sign difference between the two terms 
on the RHS of p4l) . This is necessary in order for both 
to yield an attractive force oc —V/ on a test particle in 
the shock wave background. The first term is just the 
5D shock wave solution of The second term comes 
from the residual 4D graviton zero mode which persists 
on the normal branch when g — )■ cx3 because the bulk 
volume remains finite. Its perturbative coupling is set by 
the effective 4D Planck scale = 4M|/3iJ, consis- 

tent with 1)14(1 . At all sub-horizon distances the inverse 
power of TZ in H14|l wins over the logarithm and so in 
this limit gravity really looks five-dimensional inside the 
cosmological horizon. 

The shock wave solutions derived here give us a clear 
and explicit demonstration that the gravitational "filter" 
mechanism of T] may in fact work beyond perturbation 
theory. The key is that the coefficients in ((HJ decrease 
with the index / of the Legendre polynomials. Since the 
momentum of a graviton mode on a transverse 2-sphere 
along de Sitter brane is proportional to the index of the 
polynomial, q ^ HI, this means that the modes with 
q > gH ~ 1/rc have amplitudes suppressed by q. This 
controls the rate of divergence of the series, and limits it 
to be no worse than a logarithm for any finite value of g. 
Hence at short distances TZ < rc — /g, gravitational 
shock waves are exactly the same as in 4D GR. Only for 
very low momenta are the amplitudes unsuppressed, and 
so the shock wave profile changes towards 5D at large 



distances TZ > /g. We see this explicitly in the limit 
g — > oo on the normal branch: we recover 5D gravity 
since when Tc — > all modes with finite momenta remain 
unsuppressed. Thus we expect that in DGP with finite 
g the leading order Planckian scattering behaves as in 
4D GR, since the differences are suppressed by powers 
of m/p and HTZ, both very small for highly energetic 
particles at short distances. 

We note that our exact shock waves may be a new 
arena to explore dynamics of the scalar graviton in DGP. 
Imagine probing the gravitational field of an arbitrary 
mass on the brane with a very relativistic probe. Since 
local physics in DGP obeys the usual 4D diffeomorphism 
invariance, we can transform all of the relevant physics 
to the rest frame of the probe. In this frame, the source 
mass will appear to move with a very high speed v, and 
so its gravitational field should be well approximated by 
our shock wave solutions, with the corrections due to the 
mass suppressed by the powers oi m/p = -^/l/v^ — 1. 
One can then treat the rest mass of the source as a per- 
turbation of the shock wave geometry, and study how the 
scalar graviton responds to it. If organized as an expan- 
sion in powers of m/p, perturbation theory may be under 
control. It would be interesting to investigate what hap- 
pens with the strong coupling of 0,0,0101 in this case. 
This may shed new light on the interplay of background 
curvature and scalar graviton effective field theory. 
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